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Abstract 

In this work we ofTer a significant improvement on the previous smallest spiking neural P systems 
and solve the problem of finding the smallest possible extended spiking neural P system. Paun 
and Paun [15] gave a universal spiking neural P system with 84 neurons and another that has 
extended rules with 49 neurons. Subsequently, Zhang et al. [18] reduced the number of neurons 
used to give universality to 67 for spiking neural P systems and to 41 for the extended model. 
Here we give a small universal spiking neural P system that has only 17 neurons and another 
that has extended rules with 5 neurons. All of the above mentioned spiking neural P systems 
suffer from an exponential slow down when simulating Turing machines. Using a more relaxed 
encoding technique we get a universal spiking neural P system that has extended rules with 
only 4 neurons. This latter spiking neural P system simulates 2-counter machines in linear time 
and thus suffer from a double exponential time overhead when simulating Turing machines. We 
show that extended spiking neural P systems with 3 neurons are simulated by log-space bounded 
Turing machines, and so there exists no such universal system with 3 neurons. It immediately 
follows that our 4-neuron system is the smallest possible extended spiking neural P system that 
is universal. Finally, we show that if we generalise the output technique we can give a universal 
spiking neural P system with extended rules that has only 3 neurons. This system is also the 
smallest of its kind as a universal spiking neural P system with extended rules and generalised 
output is not possible with 2 neurons. 
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1. Introduction 

Spiking neural P systems (SN P systems) [5] are quite a new computational model that are a 
synergy inspired by P systems and spiking neural networks. It has been shown that these systems are 
computationally universal [5]. Recently, Paun and Paun [15] gave two small universal SN P systems; 
They give an SN P system with 84 neurons and an extended SN P system with 49 neurons (that 
uses rules without delay). Paun and Paun conjectured that it is not possible to give a significant 
decrease in the number of neurons of their two universal systems. Zhang et al. [18] offered such a 
significant decrease in the number of neurons used to give such small universal systems. They give 
a universal SN P system with 67 neurons and another, which has extended rules (without delay), 
with 41 nciirons. Here wc give a small universal SN P system that has only 17 neurons and another, 
which has extended rules (without delay), with 5 neurons. Using a more relaxed encoding we get a 
universal SN P system that has extended rules (without delay), with 4 neurons. Table 1 gives the 
smallest universal SN P systems and their respective simulation time and space overheads. Note 
from Table 1 that, in addition to its small size, our 17- neuron system uses rules without delay. The 
other small universal SN P systems with standard rules [15,18] do not have this restriction. 

In this work we also show that extended SN P systems with 3 neurons and generalised input arc 
simulated by log-space bounded Turing machines. As a result, it is clear that there exists no such 
universal system with 3 neurons, and thus our 4-neuron system is the smallest possible universal 
extended SN P system. Following this, wc show that if we generalise the output technique we can 
give a universal SN P system with extended rules that has only 3 neurons. In addition, we show 
that a universal SN P system with extended rules and generalised output is not possible with 2 
neurons, and thus our 3-ncuron systems is the smallest of its kind. 

From a previous result [13] it is known that there exists no universal SN P system that simulates 
Turing machines in less the exponential time and space. It is a relatively straightforward matter 
to generalise this result to show that extended SN P systems suffer from the same inefficiencies. 
It immediately follows that the universal systems we present here and those found in [15,18] have 
exponential time and space requirements. However, it is possible to give a time efficient SN P 
system when we allow exhaustive use of rules. A universal extended SN P system with exhaustive 
use of rules has been given that simulates Turing machines in linear time [12]. Furthermore, this 
system has only 10 neurons. SN P systems with exhaustive use of rules were originally proved 
computationally universal by loncscu ct al. [4]. However, the technique used to prove universality 
suffered from an exponential time overhead. 

Using different forms of SN P systems, a number of time efficient (polynomial or constant time) 
solutions to NP-hard problems have been given [2,8,9]. All of these solutions to NP-hard problems 
rely on families of SN P systems. Specifically, the size of the problem instance determines the 
number of neurons in the SN P system that solves that particular instance. This is similar to 
solving problems with circuits families where each input size has a specific circuit that solves it. 
lonescu and Sburlan [6] have shown that SN P systems simulate circuits in linear time. 

In Section 4 we give a definition for SN P systems, explain their operation and give other relevant 
technical details. In Section 3 we give a definition for counter machines and we also discuss some 
notions of universality. Following this, in Section 4 we give our small universal SN P systems and 
show how their size can be reduce if we use a more relaxed encoding. In Section 5 we give our 
proof showing that extended SN P systems with 3 neurons and generalised input are simulated by 
log-space bounded Turing machines. Section 5 also contains our universal 3-neuron system with 
generalised output. We end the paper with some discussion and conclusions. 
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Section 4 
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extended! 


no 


Section 4 


4 


double-exponential 
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no 


Section 4 
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double-exponential 
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Table 1 

Small universal SN P systems. The "simulation time" column gives the overheads used by each system when simulat- 
ing a standard single tape Turing machine. ! indicates that there is a restriction of the rules as delay is not used and 
I indicates that a more generalised output technique is used. *The 18 neuron system is not explicitly given in [14]; 
it is however mentioned at the end of the paper and is easily derived from the other system presented in [14]. Also, 
its operation and its graph were presented in [11]. 

2. SN P systems 

Definition 1 (Spiking neural P system) 

A spiking neural P system (SN P system) is a tuple H = (O, ui, C72, • • • , am, syn, in, out), where: 

(i) O = {s} is the unary alphabet (s is known as a spike), 

(ii) (7i, (72, • • • , <Jm are neurons, of the form, ai = (n^, 1 ^ i ^ m, where: 

(a) Ui ^ is the initial number of spikes contained in Oi, 

(b) Ri is a finite set of rules of the following two forms: 

(i) E / s^ — > s; d, where E is a regular expression over s, b ^ 1 and d ^ 0, 

(ii) s^ X, where A is the empty word, e ^ 1, and for all E/s^ s;d from Ri 

^ L{E) where L{E) is the language defined by E, 
(Hi) syn C {1, 2, • • • , m} x {1, 2, • ■ • , m} is the set of synapses between neurons, where i ^ j for 
all {i,j) € syn, 

(iv) in, out e {u\,ai, ■ ■ ■ ,Ujn) ore the input and output neurons, respectively. 

A firing rule r = E/s^ s;d is applicable in a neuron ai if there are j ^ b spikes in a^ and 
s^ S L{E) where L{E) is the set of words defined by the rcgtilar expression E. If, at time t, rule r 
is executed then b spikes are removed from the neuron, and at time t + d the neuron fires. When 
a neuron a^ fires a spike is sent to each neuron aj for every synapse {i,j) in II. Also, the neuron 
Oi remains closed and docs not receive spikes until time t + d and no other rule may execute in a^ 
until time t + d + 1. A forgeting rule r' = s^ ^ X is applicable in a neuron a^ if there are exactly 
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e spikes in cTj. If r' is executed then e spikes are removed from the neuron. At each timestep t a 
rule must be applied in each neuron if there is one or more applicable rules at time t. Thus, while 
the application of rules in each individual neuron is sequential the neurons operate in parallel with 
each other. 

Note from 2b(i) of Definition 1 that there may be two rules of the form E/ —5- s; d, that are 
applicable in a single neuron at a given time. If this is the case then the next rule to execute is 
chosen non-deterministically. 

An extended SN P system [15] has more general rules of the form E/ s'^ — > s^; d, where 6 ^ p ^ 1. 
Thus, a synapse in an SN P system with extended rules may transmit more than one spike in a 
single timestep. The SN P systems we present in this work use rules without delay, and thus in the 
sequel we write rules as E / — > s'p . Also, if in a rule E = then we write the rule as .s'' — > .s''. 

In the same manner as in [15], spikes arc introduced into the system from the environment by 
reading in a binary sequence (or word) w e {0,1} via the input neuron ai. The sequence w is 
read from left to right one symbol at each timestep and a spike enters the input neuron on a given 
timestep iff the read symbol is 1. The output of an SN P system 11 is the time between the first 
and second firing rule applied in the output neuron and is given by the value Il{w) € N. 

A configuration c of an SN P system consists of a word w and a sequence of natural numbers 
(ri, r2, . . . , r^) where ri is the number of spikes in at and w represents the remaining input yet to 
be read into the system. A computation step cj h Cj+i is as follows: each number is updated 
depending on the number of spikes neuron ai uses up and receives during the synchronous appli- 
cation of all applicable rules in configuration Cj . In addition, if w 7^ A then the leftmost symbol of 
w is removed. A SN P system computation is a finite sequence of configurations ci , C2 , . . . . C( that 
ends in a terminal configuration Ct where for all j < t, Cj h Cj+i. A terminal configuration is a 
configuration where the input sequence has finished being read in via the input neuron (i.e. w = X 
the empty word) and either there is no applicable rule in any of the neurons or the output neuron 
has spiked exactly v times (where u is a constant independent of the input). 

Let be the x*'^ n-ary partial recursive function in a Godel enumeration of all n-ary partial 
recursive functions. The natural number value (f)x{yi,y2, ■ ■ -yn) is the result given by (f)^ on input 
(2/1,2/2, •••yn)- 

Definition 2 [Universal SN P system] A SN P system 11 is universal if there are recursive functions 
g and f such that for all x,y S N we have 4>x{yi,y2, ■■■yn) = /(n(g(x, yi, 2/2, • ■■Vn)))- 

In the next section we give some further discussion on the subject of definitions of universality. 
3. Counter machines 

Definition 3 (Counter machine) A counter machine is a tuple C = {z, R,Cm,Q,qi.qh), where 
z gives the number of counters, R is the set of input counters, Cm is the output counter, Q = 
{91,92,- •■ ,Qh} is the set of instructions, and qi,qh S Q are the initial and halt instructions, 
respectively. 

Each counter Cj stores a natural number value y 5^ 0. Each instruction qi is of one of the following 
two forms qi : INC{j),qi or q^ : DEC{j),qi,qk and is executed as follows: 
- qi : INC{j), qi increment the value y stored in counter cj by 1 and move to instruction qi. 
~ qi : DEC{j), qi, qk if the value y stored in counter Cj is greater than then decrement this value 
by 1 and move to instruction qi, otherwise if y = move to instruction g^. 
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At the beginning of a computation the first instruction executed is qi. The input to the counter 
machine is initially stored in the input counters. If the counter machine's control enters instruction 
Qh, then the computation halts at that timestep. The result of the computation is the value y stored 
in the output counter when the computation halts. 

We now consider some different notions of universality. Korec [7] gives universality definitions 
that describe some counter machines as weakly universal and other counter machines as strongly 
universal. 

Definition 4 [Korec [7]] A register machine M will be called strongly universal if there is a recur- 
sive function g such that for all a;,y G N we have ^x{y) = ^1^(5(2;), y)- 

Here $^((/(a;), y) is the value stored in the output counter at the end of a computation when M is 
started with the values g{x) and y in its input counters. Korec's definition insists that the value y 
should not be changed before passing it as input to M. However, if we consider computing an n-arry 
function with a Korec-strong universal counter machine then it is clear that n arguments must be 
encoded as a single input y. Many Korec-strong universal coimter machines would not satisfy a 
definition where the function (f>x in Definition 4 is replaced with an n-arry function with n > 1. 
For example, let us give a new definition where we replace the equation ''(j)x{y) — ^^/(^(x), y)" 
with the equation "^^(j/i, 2/2, • • • , J/n) = ^^^ig{x),yi,y2, ■ ■ ■ ,yn)" in Definition 4. Note that for 
any counter machine M with r counters, if r < n then M does not satisfy this new definition. 
It could be considered that Korec's notion of strong universality is somewhat arbitrary for the 
following reason: Korec's definition will admit machines that require n-arry input (yi, 2/2, • • • , 2/n) 
to be encoded as the single input y when simulating an n-arry function, but his definition will not 
admit a machine that applies an encoding function to y (e.g. y^ is not permitted). Perhaps when 
one uses this notion of universality it would be more appropriate to refer to it as strongly universal 
for unary partial recursive functions instead of simply strongly universal. 

Korec [7] also gives a mimbcr of other definitions of universality. If the equation (pxiy) = 
^\f{g{x),y) in Definition 4 above is replaced with any one of the equations <j)x{y) = ^M{92{x,y)), 
(pxiv) = /(^M(ff(^)' y)) 4'x{y) = f{^M{92{x,y))) then the counter machine M is weakly uni- 
versal. Korec gives another definition where the equation (j)x{y) = ^\i{g(x),y) in Definition 4 is 
replaced with the equation 0x(y) = /(^m(5(^)' ^iv)))- However, he does not include this definition 
in his list of weakly universal machines even though the equation (j)x{y) = f{'^\i{9{x),h{y))) allows 
for a more relaxed encoding than the equation (j)x{y) = f{^%[{g{x),y)) and thus gives a weaker 
form of universality. 

For each number m > 2 there exists universal m-counter machines that allow </>" and its input 
(2/1, 2/2, • • • , J/n) to be encoded separately (e.g. via g{x) and /i"(yi, y2, ?/«))• For universal 2- 
counter machines all of the current algorithms encode the function and its input (?/i, ?/2, ■ ■ ■ , Vn) 
together as a single input (e.g. via 9"+^(.'i:, r/i, y2, ■ ■ ■ , yn))- Using such encodings it is only possible to 
give universal 2-counter machines that Korec would class as weakly universal. Some other limitations 
of 2-counter machines were shown independently by Schroeppel [16] and Barzdin [1]. In both cases 
the authors are examining unary functions that are uncomputable for 2-counter machines when the 
input value; to the counter machine must equal the input to the function. For example Schroeppel 
shows that given n as input a 2-counter machine cannot compute 2". It is interesting to note that 
one can give a Korec-strong universal counter machine that is as time/space inefficient as a Korec- 
weak universal 2-countcr machine. Korec's definition of strong universality deals with input and 
output only and is not concerned with the (time/space) efficiency of the computation. 
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In earlier work [15], Korec's notion of strong universality was adopted for SN P systems^ as 
follows: A spiking neural P system 11 is strongly universal if n(10^^^10^~^l) ~ (f>x{y) for all x and 
y (here if (j)x{y) is undefined so to is n(10^~^10'^~^l)). As with the SN P systems given in [15,18], 
the systems we give in Theorems 1 and 3 satisfy the notion of strong universality adopted from 
Korec in [15]. Analogously, our system in Theorem 2 could be compared to what Korec refers to as 
weak universality. However, as we noted in our analysis above, it could be considered that Korec's 
notion of strong universality is somewhat arbitrary and we also pointed out some inconsistency in 
his notion of weak universality. Hence, in this work we rely on time/space complexity analysis to 
compare the encodings used by small SN P system (see Table 1). 

It is well known that counter machines require an exponential time overhead to simulate Turing 
machines [3]. Counter machines with only 2 counters arc \mivcrsal [10], however, they simulate 
Turing machines with a double exponential time overhead. In the seqiiel we give some universal 
SN P systems that simulate 3-counter machines and others that simulate 2-counter machines. The 
reason for this is that when using our algorithm there is a trade-off between the size and the 
time efficiency of the system. This trade-off is dependant on whither we choose to simulate 3- 
counter machines or 2-counter machines. When simulating Turing machines, 3-counter machines 
suffer from an exponential time overhead and 2-counter machines suffer from a double-exponential 
time overhead, and thus the simulation of 3-counter machines is preferable when considering the 
time efficiency of the system. If we are considering the size of our system then 2-countcr machines 
have an advantage over 3-counter machines as our algorithms require a constant number of neurons 
to simulate each counter. 

4. Small universal SN P systems 

We begin this section by giving our two extended universal systems Ylc^ and Hc^ , and following 
this we give our standard system H^^ . We prove the universality of Hcs and H^^ by showing that 
they each simulate a universal 3-counter machine. From Hcg we obtain the system U'q^ which 
simulates a universal 2-counter machine. 

Theorem 1 Let C3 be a universal counter machine with 3 counters that completes it computation 
in time t to give the output value Xo when given the pair of input values (xi, X2)- Then there is a 
universal extended SN P system Hcg that simulates the computation of C3 in time 0{t+xi +X2+X0) 
and has only 5 neurons. 

PROOF. Let C3 = (3, {ci, C2}, C2, Q, qi,qh) where Q = {qi,q2, ■ • ■ , qh}- Our SN P system Hcg is 
given by Figure 1 and Table 4. The algorithm given for Hc^ is deterministic. 

4.0.1. Encoding of a configuration of C:^ and reading input into H^g 

A configuration of C3 is stored as spikes in the neurons of Hcg . The next instruction qi to be 
executed is stored in each of the neurons 0-2, 0-3 and 0-4 as 4(ft, + i) spikes. Let x\, X2 and a;3 be the 
values stored in counters ci, C2 and C3, respectively. Then the values xi, X2 and x^ are stored as 
&h{xi -\- 1), 8h{x2 + 1) and 8h{x3 + 1) spikes in neurons (72, as and 0-4, respectively. 

The input to Hcg is read into the system via the input neuron cti (see Figure 1). If C3 begins its 
computation with the values xi and X2 in coimters ci and C2, respectively, then the binary sequence 
w = 10^^~^10^^^^1 is read in via the input neuron o"i. Thus, cti receives a single spike from the 

^ Note that no formal definition of this notion was explicitly given in[15]. 
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Fig. 1. Universal extended SN P system Tlc^- Eadci oval labeled (7i is a neuron. An arrow going from neuron u, to 
neuron illustrates a synapse 

environment at times ti, ^xi+i and txi+x2+i- We explain how the system is initiahsed to encode 
an initial configuration of C3 by giving the number of spikes in each neuron and the rule that is 
to be applied in each neuron at time t. Before the computation begins neuron ai initially contain 
8h spikes, as contains 2 spikes, 174 contains 8h + 1 spikes and all other neurons contain no spikes. 
Thus, when cti receives it first spike at time ti we have 

h : 

ai= 8/1 + 1, s^^+'/s^f'^s^^, 

= 2, s^/s — >■ s, 

a4 = 8/1+1, s^^+V/" ^ .s^"-!. 

where on the left dk = z gives the number z of spikes in neuron crfe at time t and on the right is the 
rule that is to be applied at time t, if there is an applicable rule at that time. Thus, from Figure 1, 
when we apply the rule s^^j+i^gSft, _^ ^sh -^^ neuron ai, s^/s ->• s in 0-3, and gS/i+i^gSft, _^ gS'^-i in 
(74 at time ti we get 

o\ 

0-2 
0-3 
0-4 
0-5 

ts : 

ci 
01 
0-3 

(74 
0-5 



8/1+1, s8'>+l/s8/> ^ //.^ 
8/1, 

8/1 + 1, s«''+Vs^'* ^ s, 

8/1 + 1, s8/»+l/g8/. ^ ^8?»-l^ 

1, s ^ A, 



8/1 + 1, s«''+Vs^'* ^ 

16/1, 

8/1+1, s8''+Vs«'' ^ s, 

8/1+1, s8?.+l/g8/. ^ ^8?.-l^ 

1, s -)■ A. 
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Neuron ai fires on every timestep between times ti and t^i+i to send a total of Shxi spikes to a2 
thus we get 



t 



Xl + l 



ai=8h + 2, 
CT2 = 8hxi, 

£73 = 8/l + 1, 

(74 = 8/l + 1, 
0-5 = 1, 



Sh+2 Ig&h+l _^ ^S,h+1 



s A, 



txi+2 ■ 

(Ti = 8/l+ 1, 

(72 = 8h{xi + 1) + 1, 
(73 = 8/i + 2, 
(74 = 8/i + 2, 

0-5 = 1, 



Rh + 1 /„8ft 8/l 
r 5 1 



,8/1+2 /^8/i ^ c8/j-l 

s A, 



txi+3 ■ 

ai = 8ft + 1, 

(72 = 8/i(a;i + 1) + 1, 

£73 = 16h + 2, 

(74 = 8/l + 2, 



Neuron cri fires on every timestep between times t^^^i and txi+x2+i to send a total of 8hx2 spikes 
to (73. Tlius, wlien ai receives tlie last spike from its environment we have 



ai=8h + 2, 

(72 = 8h{xi + 1) + 1, 
(73 = 8hx2 + 2, 
(74 = 8/i + 2, 



^8/»+2 , 8/. ^ ^8fc-l 



*xi+x2+2 : 

(71 = 8/l + 1, 

(72 = 8/i(a;i + 1) + 2, 
(73 = 8h{x2 + 1) + 3, 

(74 = 8/l + 3, 



^^8/.)*^8/.+2/^8/.+2 ^ ^2/.^ 
^8h+3 ^ ^2/._ 
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ixi+X2+3 • 

(Ti = 6h + 1, 
0-2 = 8h{xi + 1), 
0-3 = Sh{x2 + 1), 
(74 = 8h, 
(75 = 2/l, 



s^'' ^ A, 



txi+X2+4 '■ 

a2 = 8/i(xi + l) + 4(/i+l), 
(73 = 8/i(a;2 + l)+4(/i+l), 
0-4 = 8/i + 4(/i+ 1). 

At time txi+x2+4: neuron (T2 contains 8ft-(a:;i + 1) + 4(ft, + 1) spikes, as contains 8/1(2:2 + 1) + 4(ft, + 1) 
spikes and a4 contains 8h + 4(/i + 1) spikes. Thus at time txi+x2+4 the SN P system encodes an 
initial configuration of C3. 



4.0.2. Ilcg simulating qi : INC{\),qi 

Let counters Ci, C2, and C3 have values Xi, X2, and 0:3, respectively. Then the simulation of 
qi : INC{l),qi begins at time tj with 8/i(xi + 1) + 4(/i + i) spikes in (72, 8h{x2 + 1) + 4:{h + i) spikes 
in 1J3 and 8h{x3 + 1) + 4(/i + i) spikes in 0-4. Thus, at time tj we have 



(72 = 8/i(xi + 1) + 4(/i + i), (s8'»)*s4('i+») /s*(h+i) 



s 



4(h+i) 



(73 = 8h{x2 + 1) + 4(/l + i), (s8'')*s4(''+*) /g8fc+4(/,+i) ^ g6/»^ 

(74 = 8h{xs + 1) + 4(/l + i), (^gShyg4(h+i) l^Sh+iih+i) ^ g6/._ 

Prom Figure 1, when we apply the rule (^g^hy gi{h+i) j gi{h+i) _^ gi(h+i) neuron (72 and the rule 

(^gihygiih+i) igih+A{h+i) _^ g6h tijjjg ^. 

ai = 16/i + 4z, gi6h+4i ^ ^m+4;^ 
(72 = 8h{xi + 1), 

(73 = 8/1X2, 
(74 = 8/1X3, 

(75 = 6/1, s'^'' A, 



tj+2 : 

(72 = 8/i(a;i+2) + 4(/i + 0, 

(73 = 8h{x2 + 1) +4(/i + /), 
(74 = 8h{x3 + l)+4{h + l), 



At time tj+2 the simulation of (/j : INC{l),qi is complete. Note that an increment on the value xi 
in counter ci was simulated by increasing the 8h{xi + 1) spikes in (72 to 8h{xi + 2) spikes. Note also 
that the encoding 4(/i + 1) of the next instruction qi has been established in neurons (72, as and a^. 
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4.0.3. 11(73 simulating qt : DEC{l),qi,qk 

There are two cases to consider here. Case 1: if counter ci has value Xi > 0, then decrement 
counter 1 and move to instruction g^+i. Case 2: if counter ci has value xi = 0, then move to 
instruction qk- As with the previous example, our simulation begins at time tj. Thus Case 1 (a;i > 
0) gives 



(72 = 8h{xi + 1) +4(/i + i), 
0-3 = 8h{x2 + 1) + 4(/i + i), 
0-4 = 8h{x3 + 1) + 4(/i + i), 



.Shy gl6h+4{h+i) ^ gl2h+4i _^ g6h+4i 
{S^hyg4ih+i)/^4ih+i) ^g2h^ 
f^gShy g4:(h+i) lgi(h+i) _^ g2h^ 



tj+i : 

(71 = lOh + 4i, 
(72 = 8hxi, 

(73 = 8hix2 + 1), 
(74 = 8h{x3 + 1), 
(75 = 2h, 



glOh+4i _^ g4ih+l) 



tj+2 ■■ 

(72 = 8hxi +4(/i + I), 

(73 =8/l(x2 + l)+4(/l + 0, 

a4 = 8/i(a;3 + l)+4(/i + 0. 



At time tj-^-2 the simulation of qi : DEC{1), qi, qk for Case 1 (xi > 0) is complete. Note that a 
decrement on the value xi in counter ci was simulated by decreasing the 8h{xi + 1) spikes in (72 to 
8hxi spikes. Note also that the encoding 4(/i + /) of the next instruction qi has been established in 
neurons (72, (73 and (74. Alternatively, if we have Case 2 = 0) then we get 



0-2=8/1 + 4(/i + i), 

(73 = 8/1(0:2 + 1) + 4(/i + i), 

(74 = 8/i(a;3 + 1) + 4(/i + i). 



g8/i+4(/(+i) Ig4(h+i) _^ g4(h+i) ^ 
(^gShy g4(h+i) Ig4(h+i) _^ g2h^ 
f^gShy g4{h+i) Ig4(h+i) _^ g2h^ 



(7i = 8/l + 4i, 
(72 = 8h, 
(73 = 8/i(a;2 + 1), 
(74 = 8/l(X3 + 1), 

(75 = 2/l, 



^8/1+41 _^ g4{h+k) 
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tj+2 : 

(72 = 8/i + 4(/i + fc), 

(73 = 8h{x2 + 1) +4(/i + fc), 

(T4 = 8/i(a;3 + 1) + 4(/i + fc). 

At time the simulation of : DEC{1), qi, qk for Case 1 (xi = 0) is complete. The encoding 
A[h + k) of the next instruction q^ has been established in neurons (72, cts and (74. 



4.0.4. Halting 

The halt instruction g^j is encoded as 4/i + 5 spikes. Thus, if C3 enters the halt instruction we 
get 

(72 = 8/i(a;i + 1) + 4/i + 5, 

£73 = 8/1(3:0 + 1) + 4/l + 5, (g8/»^*g20fc+5/gl2h ^ 

(74 = 8/i(a;3 + 1) + 4/i + 5, 
CTi =2, 

(72 = 8/i(a;i + 1) +4/1 + 5, 

(73 = 8/1x0 + 5, 
(74 = 8/i(x3 + 1) + 4/1+ 5, 

= 2, 



ij+2 : 

(7i = 1, s ^ A, 

(72 = 8/i(xi + 1) + 4/i + 5, 

(73 = 8/i(xo - 1) + 5, (s«'')*s'®''+Vs*'' ^ s, 

(74 = 8/l(X3 + 1) + 4/1 + 5, 

cr5 = 1, S — >■ A. 

The rule {^g^hy gWh+5 j g8h _^ § jg applied a further Xo — 2 times in (73 until we get 

(71 = 1, s — > A, 
0-2 = 8/i(a;i + 1) + 4/1 + 5, 

(73 = 8/1 + 5, s^''+^^•s^ 
(74 = 8/1(0:3 + 1) + 4/1 + 5, 

(75 = 1, s — >^ A. 
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tj+x„ + l ■ 

ai =2, -)■ A, 

0-2 = 8h{xi + 1) + 4/1 + 5, 
tT4 = 8/;,(a;3 + 1) + 4/i + 5, 

(75 — 2, — >• S. 

As usual the output is the time interval between the first and second spikes that are sent out of the 
output neuron. Note from above that the output neuron fires for the first time at timestep tj^i 
and for the second time at timestep tj^xo+i - Thus, the output of Ilca is Xo the value of the output 
counter C2 when C3 enters the halt instruction qh- Note that if a;2 = then the rule s^^^~^^ — )• is 
executed at timestep tj, and thus only one spike will be sent out of the output neuron. 

We have now shown how to simulate arbitrary instructions of the form qi : INC{l),qi and 
qi : DEC{l),qi,qk that operate on counter ci. Instructions which operate on counters C2 and C3 
are simulated in a similar manner. Immediately following the simulation of an instruction Ilc^ is 
configured to simulate the next instruction. Each instruction of C3 is simulated in 2 timesteps. The 
pair of input values (xi , X2 ) is read into the system in Xi + a;2 + 4 timesteps and sending the output 
value Xo out of the system takes Xo + 1 timesteps. Thus, if C3 completes it computation in time t, 
then simulates the computation of C3 in linear time 0{t + xi + X2 + Xo)- □ 

Theorem 2 Let C2 he a universal counter machine with 2 counters thai completes it computation 
in time t to give the output value Xo when given the input value xi. Then there is a universal 
extended SN P system XIcj that simulates the computation of C2 in time 0{t + xi + Xo) and has 
only 4 neurons. 

PROOF. Let C2 = (2,{ci},C2,Q,9i,g/i) where Q = {qi,q2,--- ,Qh}- The rules for the SN P 

system Tlc2 a-rc given by Table 5 and a diagram of the system is obtained by removing neuron 
(74 from Figure 1. If C2 begins its computation with the value xi in counter ci then the binary 
sequence w = 10^^~^1 is read in via the input neuron ai. Before the computation begins neurons 
(7i, (72, (73 and (J5 respectively contain 8h, 8h + 1, 16h + 1 and spikes. Like Tlc^, ^0-2 encodes the 
value X of each counter as 8h{x + 1) spikes and encodes each instruction qi as 4(/i + i) spikes. The 
operation of IIcj is very similar to the operation of Hc^ , and thus it would be tedious and repetitive 
to go through another simulation here. Tlc2 simulates a single instruction of C2 in 2 timesteps in a 
manner similar to that of flcg- The inputting and outputting techniques, used by IIcj, also remain 
similar to those of Ilcg, and thus the running time of is 0{t + Xi + Xo). □ 

The SN P system in Theorem 3 simulates a counter machine with the following restriction: if a 
counter is being decremented no other counter has value at that timestep. Note that this does 
not result in a loss of generality as for each standard counter machine there is a counter machine 
with this restriction that simulates it in linear time without an increase in the number of counters. 
Let C be any counter machine with m counters. Then there is a counter machine C with m 
counters that simulates C in linear time, such that if C" is decrementing a counter no other counter 
has value at that timestep. Each counter in C that has value y is simulated by a counter in 
C that has value y + 1- The instruction set of C is the same as the instruction set of C with 
the following exception each qi ; DEC{j), qi, qk instruction in C is replaced with the instructions 
: DEC{3)q'^,q[), {q[ : DEC{j)ql ,ql), (g* : INC{j),qi), and (g* : INC{j),qk). The reason we 
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Fig. 2. Part 1 of the universal SN P system H^^. Each oval labeled ai is a neuron. An arrow going from neuron <7i 
to neuron dj illustrates a synapse 

need these extra instructions is that y is encoded as y + 1 and we must decrement twice if we wish 
to test for an encoded 0. 

Theorem 3 Let C3 be a universal counter machine with 3 counters and h instructions that com- 
pletes it computation in time t to give the output value Xg when given the input (.xi,X2). Then there 
is a universal SN P system Wq^ that simulates the computation of C3 in time 0{ht + Xi + X2 + Xg) 
and has only 17 neurons. 

PROOF. Let C3 = {3,{ci,C2},C3,Q,qi,qh) where Q = {qi,q2,--- ,qh}- Also, without loss of 
generality wc assume that during Cs's computation if C3 is decrementing a counter no other counter 
has value at that timestep (sec the paragraph before Theorem 3). The SN P system 11^^ is given 
by Figures 2 and 3 and Tables 6 and 7. As a complement to the figures, Table 3 may be used to 
identify all the synapses in 11^^ . The algorithm given for 11^^ is deterministic. 

4.0.5. Encoding of a configuration of Cz and reading input into 11^^ 

A configuration of C3 is stored as spikes in the neurons of H^^. The next instruction qi to be 
executed is stored in each of the neurons £72, (T3, (T4, 0-5, erg) and (T7 as 21{h + i) + 1 spikes. Let xi, 
X2 and X3 be the values stored in counters ci, C2 and C3, respectively. Then the value xi is stored 
as 6(a;i + 1) spikes in neuron as, X2 is stored as 6(2:2 + 1) spikes in ag, and 0:3 is stored as 6(0:3 + 1) 
spikes in aio- 
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The input to IIq^ is read into the system via the input neuron cti (see Figure 2). If C3 begins 
its computation with the values xi and X2 in counters ci and C2, respectively, then the binary 
sequence w = 10^'^~^10^^~^1 is read in via the input neuron cti. Thus, cti receives a spike from the 
environment at times ti, tx^+i and txi+x2+i- We explain how the system is initialised to encode an 
initial configuration of C3 by giving the number of spikes in each neuron and the rule that is to be 
applied in each neuron at time t. Before the computation begins neurons fT2, tra, (74, (T5, ae and (77 
each contain 40 spikes, neurons erg, erg and aio each contain 3 spikes, and neurons (T12, (T13 and 0-14 
each contain 21h — 2 spikes. Thus, when ai receives it first spike at time ti we have 

h : 

(Tl = 1, s — >■ s, 

O'2,Cr3,Cr4,CT5,cr6>O'7 = 40, 

as,ag,aio = 3, 

cTi2,ai3,ai4 = 21/i-2, {s^Ys^/s^^s. 

Thus, from Figures 2 and 3, when we apply the rule s s in neuron cti and the rule {s^)*s^/s^ — >■ s 
in 0-12, and 0-14 at time ti we get 

t2 : 

O'2,CT3,O-4,0-5,a-6,a'7 = 41, 

as,(T9,aio = 4, 
0-12, 0-13,0-14 = 21h - 4, 
0-15, 0-16, 0-17 = 3, 

is : 

0-2, 0-3, 0-4, 0-5, 0-6, 0-7 = 41, 

<J8 = 10, 

09, 0-10 = 10, 
0-11 = 6, 

0-12,0-13,0-14 = 21/1 - 4, 
0-15, 0-16, Oi7 = 3. 

ti : 

0-2, 0-3, 0-4, 0-5, 06, 0-7 = 43, 

(78 = 16, 

0-9, oio = 10, 
0-11 = 7, 

0-12,013,014 = 21/1 - 4, 

0-15, 0-16, Oi7 = 3. 

Neurons 02, 0-3, 0-4, 05, erg and 07 fire on every timestep between times t2 and txi+2 to send a total 
of 6x1 spikes to os, and thus we get 



s /s -)■ s, 



s /s -)■ s, 
^ A, 



A, 
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txi+1 ■ 

(Jl = 1, S — >■ S, 

0-2, 0-3, 0-4, 0-5, CT6, ay = 43, s'^^/s^^s, 
CT8 = 6(a;i - 1) + 4, 

a9,aio = 10, (/)*si7s6^s, 

(Til =7, -s- A, 

0-12, 0-13,(714 = 21h - 4, 

t2:i+2 : 

(72 = 44, s^Vs^^ ^ s> 

0-3, 174, 0-5, 0-6, (77 = 44, S^^/s^^-)-S, 

(78 = 6x1+5, (s*^)*s"/s*^ ^ s, 

(7g = 11, 

(710 = 11, (sS)*s"/s' 

(711 =7, ^ A, 

0'12,(7l35Cri4 = 21/1 — 3, 
Cl5,(7l6jCri7 = 3, 

ia:i+3 : 

(72 = 22, S^^/S^ -)■ S, 

(73, £74, £75, £76, (77 = 16, S^^/s^->-S, 

(78 = 6X1 + 5, (S^)*s"/S^ S, 

erg = 17, 

(710 = 11, (s'^)*s"/s' ^ s, 

(711 =7, ^- A, 

0'12, 1713, £714 = 21/1 — 3, 
Cl5jCri6,£7l7 = 3. 



Neurons U2, (73, (74, CT5, CTg and a-j fire on every timestep between times txi+2 and tx-i+x2+2 to send 
a total of 6x2 spikes to £79. Thus, when £7i receives the last spike from its environment we have 
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ai = 1, s — >■ s, 

a2 = 22, s^Vs^ s, 

(73, (74, (T5, (Tg, (77 = 16, S^^/s^^S, 

CT8 = 6xi+5, (/)*s"/s^^s, 

(79 = 6X2 + 5, 

aio = 11, {S^S^'/S^ ^ S, 

(711 =7, ^- A, 

^12, CIS, (7i4 = 21h - 3, 
^15 5 CTie, (7i7 = 3, 

txi+X2+2 ■ 

(72 = 23, 5^3/^5 ^ 

0'3,cr4,0"5,0"6,0"7 = 17, 
(78 = 6(xi + 1), 
(79 = 6{X2 + 2), 

aio = 12, 

(7ii =7, s'^ — > A, 

0-12,0-13,(714 = 21/1-2, (s3)*s4/s3 ^ s, 

0-15, 0-16, (7l7 = 3, 

txi+X2+3 '■ 

0-2, 0-3, (74, (75, (76, (77 = 18, 

(78 = 6(a;i + 1) + 1, (s'^)*siVs ^ s> 

(79 = 6(X2 + 2) + 1, (/)*s7s^ ^ s, 

(710 = 13, (/)*s7s^ ^ s, 

(711 = 1, s -s- A, 

Cri2, 0-13, 0-14 = 21/1-5, (s^)*S^/s3 S, 

0-15, 0-16, 0-17 = 6, 

txi+X2+4 '■ 

0-2,0-3,(74,(75,(76,(77 = 21, 

as = Q{xi + 1), 
0-9 = 6(X2 + 1), 
0-10 = 6, 

(711 =1, s ^ A, 

0-12,0-13,0-14 = 21/1-8, {s^)*s^/s^ ^ s, 

0-15,0-16,0-17 = 9. 
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After a ftirther 7/i — 3 timestep we get 



txi+X2+7h+l '■ 

0'2,0'3,CT4,O'5,O'6,Cr7 = 21, 

as = 6(.Ti + 1), 

(Tg = 6{X2 + 1), 
CTlO = 6, 

ai2 = 1, s ^ s, 

Cl3, cri4 = 1, S ->■ A, 

O'lS, cri6i CiT = 21ft,, 

txi+X2+7h+2 ■ 

0'2,cr3,cr4,(T5,cr6,iJ7 = 21, 

as = 6(xi + 1), 
erg = 6(a;2 + 1), 
CTlO = 6, 

0-15, 0-16, (717 = 21/1 + 1, {S^)*S^/S^ ^ S, 

txi+X2+7h+3 '■ 

a2,a3,a4,a5,ae,a7 = 21 + 3, 
as = Q{xi + 1), 

(Tg = 6{X2 + 1), 

cTio = 6, 

CTl2, CIS, (7i4 = 3, 

(715,0-16,(717 = 21/1-2, (s^)*sVs^ S. 



Neurons (715, aie and £717 continue to fire at each timestep. Thus, after a further 7h — 1 steps we 
get 



txi+X2 + 14h+2 '■ 

(72, (73, (74, (75, ae, (77 = 21/1 + 21, 
as = 6{xi + 1), 
ag = 6{X2 + 1), 
(710 = 6, 

(7l2, (713, (7l4 = 21/1, 

(7l5 = 1, S ^ S, 

(716,(717 = 1, S ^ A. 
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Fig. 3. Part 2 of the universal SN P system 11^^ . Each oval labeled (Tj is a neuron. An arrow going from neuron Oi 
to neuron Cj illustrates a synapse 

ia;i+a:2 + 14?i+3 '■ 

a2,cF3,a4,a5,a6,a7 = 21{h+ 1) + 1, 
as = G{xi + 1), 

(Tg = 6(X2 + 1), 
<7W = 6, 

0-12, 0-13,0-14 = 21h+ 1. 

At time txi+x2+i4h+3 neurons (T2, cts, 0-4, 0-5, ere and 0-7 each contain 21(/i + l) + l spikes, cts contains 
6(a;i + 1) spikes, ug contains 6{x2 + 1) spikes and ctio contains 6 spikes. Thus, at time txi+x2+i4h+3 
the SN P system encodes an initial configuration of C3. 

4.0.6. Algorithm overview 

Here we give a high level overview of the simulation algorithm used by E^^. Neurons erg, o-g and 
aio simulate the counters of Ci, C2 and C3, respectively. Neurons CT2, 0-3, 0-4. (T5, (Tq and cry are the 
control neurons. They determine which instruction is to be simulated next by sending signals to the 
neurons that simulate the counters of C3 directing them to simulate an increment or decrement. 
There are four different signals that the control neurons send to the simulated counters. Each of 
these signals takes the form of a unique number of spikes. If 1 spike is sent to erg, o-g and ctio then 
the value in ug (counter ci) is tested and ug (counter C2) and aio (counter C3) are decremented. If 
2 spikes are sent the value of (T9 is tested and as and aiQ arc decremented. If 3 spikes arc sent the 
value of aio is tested and as and erg are decremented. Finally, if 6 spikes are sent all three counters 
are incremented. Unfortunately, all of the above signals have the effect of changing the value of 
more than one simulated counter at a time. We can, however, obtain the desired result by using 
more than one signal for each simulated timestep. If we wish to simulate INC we send 2 signals 
and if we wish to simulate DEC we send either 8 or 2 signals. Table 2 gives the sequence of spikes 
(signals) to be sent in order to simulate each counter machine instru(;tion. To explain how to use 
Table 2 we will take the example of simulating INC {2). In the first timestep, all three simulated 
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IllsllUttitJll 



StXlUCllCC of S])ikcb -St!lll fl Ulll CT2; CT:; . (T4. tT5, tT(, Jlllcl (77 



/Arc(i) 
INC(2) 
INC{3) 
DEC{1) 
DEC{1) 
DEC{2) 
DEC{2) 
DEC{3) 
DEC{3) 



3, 0, 6 if X3 = 



6, 1 



6, 2 



6, 3 



3, 0, 6, 6, 6, 2, 2, 1, 1 if 13 > 



2, 0, 6 if a;2 = 



2, 0, 6, 6, 6, 3, 3, 1, 1 if a;2 > 



1, 0, 6, 6, 6, 3, 3, 2, 2 if XI > 



1, 0, 6 if xi = 



Table 2 



This table gives a counter machine instruction in the left column followed, in the right column, by the sequence that 
is used by 11^^ to simulated that instruction. Each number in the sequence represents the total number of spikes to 
be sent from the set of neurons 1T2, fa, 1T4, <75, ere and err at each timestep. 

counters erg, erg and (Tiq are incremented by sending 6 spikes, and then in the second timestep the 
simulated counters erg and aio are decremented by sending 2 spikes. This has the effect of simulating 
an increment in counter C2 and leaving the other two simulated counters unchanged. 

Each counter machine instruction qi is encoded as 21{h + i) + 1 spikes in each of the control 
neurons. At the end of each simulated timestep the number of spikes in the control neurons must 
be updated to encode the next instruction qk. The update rule s2i(/i+t)-2ife _^ ^ jg applied in 
each control neuron leaving a total of 21fc spikes in each control neuron. Following this, 21h + 1 
spikes are sent from neurons ctis, aie and air to each of the control neurons. This gives a total of 
21{h + k) + 1 spikes in each control neuron. Thus encoding the next instruction q^- (Note that the 
rule s2i(ft,+i)-2ife _^ s jg simplification of the actual rule used.) 

4.0.7. U'^^ simulating qi : INC{l),qi 

The simulation of INC{1) is given by the neurons in Figures 2 and 3. Let xi, X2 and X3 be the 
values in counters ci, C2 and C3 respectively. Then our simulation of qi : INC{l),qi begins with 
6(a;i + 1) spikes in as, 6{x2 + 1) spikes in erg, 6(a;3 + 1) spikes in aio, 21(/i + i) + 1 spikes in each of 
the neurons (72, as, a^, a^, ag and ar, and 21h + 1 spikes in each of the neurons 0-12, cria and ai4. 
Beginning our simulation at time tj, we have 



(72 = 21{h + i) + 1, 
(73,(74,(75,(76,(77 = 21(/l + i) + 1 




as = 6(xi + 1), 
(79 = 6{X2 + 1), 
o-io = 6(a;3 + 1), 



0'12,0-i3,C7i4 = 21h+ 1, 
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Thus, from Figures 2 and 3 we get 
tj+i : 



(72 = 21(/i + i) - 2, 

(73,(T4,(T5,(T6,0'7 = 211 — 4, 

as = 6(a;i + 2), 

erg = 6{X2 +2), 

(710 = 6(a;3 + 2), 
an = 6, 

0'l2,CTi3,Cri4 = 21h - 2, 
CIS, cri6, cri7 = 3, 




/ A, 



(^2^ a4, a^, a(i, a^ = 21Z — 3, 



(78 = 6(.Tl + 2) + 1, 
(79 = 6(X2 + 2) + 1, 

o-io = 6(x3 + 2) + 1, 

CTll = 1, 



s ->■ A, 



o'i2, CIS, = 21/i — 5, 



fi5, cie, ciy = 6, 

f2,0'3,C4,0'5jf6,cr7 = 21Z, 

(78 = 6(xi + 2), 
(79 = 6(a;2 + 1), 
(710 = 6(X3 + 1), 

CTll =1, s ^ A, 

0-12, (713, (714 = 21/1 -8, (s^)*sVs^ ^- S, 

fi5, cie, ciy = 9. 



The remainder of this simulation is similar to the computation carried out at the end of the 
initialisation process (see the last paragraph of Section 4.0.6 and timesteps txi+x2+4 to fxi+x2+i4/i+3 
of the Section 4.0.5). Thus, after a further 14/i — 1 timesteps we get 



tj+14:h+2 '■ 

a2,a3,a4,,a5,ae,a7 = 21{h + l) + 1, 
as = &{xi + 2), 
ag = 6{X2 + 1), 
(710 = 6(a;3 + 1), 

0-12,(713,(714 = 21/1+ 1, {s^)*s'^/s^ ^ s. 
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At time the simulation of qi : INC{1), qi is complete. Note that an increment on the value 

xi in counter Ci is simulated by increasing the number of spikes in as from 6{xi + 1) to 6(a;i + 2). 
Note also that the encoding of the next instruction qi is given by the 21{h + l) + l spikes in neurons 
o'2, cTs, fT4, 0-5, (76 and cry. 

4.0.8. H'q^ simulating qi : DEC{l),qi,qk 
If we are simulating DEC{1) then we get 

tj-- 

(72 = 21(/l + i) + l, s21(''+') + Vs^ ^ S, 

c^s, 0-4, 0-5, ere, 0-7 = 21(/i + i) + 1, 

as = 6(a;i + 1), 
erg = 6(2:2 + 1), 
(^10 = 6(a;3 + 1), 

0-12, 0-13, <Ti4 = 21/1+ 1, (S^)*S^/S^ S. 

To help simplify configurations we will not include neurons 1712, (T13, and au until the end of the 
example. When simulating DEC{1) there are two cases to consider. Case 1: if counter ci has value 
xi > 0, then decrement counter 1 and move to instruction qi+i. Case 2: if counter ci has value 
Xi = 0, then move to instruction g^. In configuration tj^i our system determines if the value xi 
in counter 1 is > by checking if the number of spikes in erg is > 13. Note that if we have Case 1 
then the rule {s^)*s^^/s — ?• s is applied in as sending an extra spike to neurons CT2, 173, (74, 0-5, ag 
and (77 thus recording that xi > 0. Case 1 proceeds as follows: 

tj+i : 

(72 = 21(/i + i) -4, 
o'3, 0-4, 0-5, ere, 0-7 = 21{h + i) + 2, 
as = 6{xi + 1) + 1, 
£79 = 6{X2 + 1) + 1, 
(710 = 6(a;3 + 1) + 1, 
CTll = 1, 

tj+2 ■ 

(72 = 21{h + i) - 1, 

(T'S, 0-4, £75, (76, (77 = 21(/l + i) + 5, 

as = 6{xi + 1), 

(79 = 6X2, 

aio = 6x3 

ail = 1, s — > A. 

The method we use to test the value of as (simulated counter ci) has the side-effect of decrementing 
(79 (simulated counter C2) and aio (simulated counter C2). Following this, in order to get the correct 
values our algorithm takes the following steps: Each of our simulated counters {as, ag and crio) are 



(s6)*siVs ^ S, 
(s<')*sVs^ ^ s, 
(s<')*sVs^ ^ s, 
s ^ A, 
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incremented 3 times, and then the simulated counter erg is decremented 4 times, whilst the simulated 
counters erg and (Jiq are each decremented twice. Thus, the overall result is that a decrement of ci 
is simulated in as and the other encoded counter values in erg and aio remain the same. Continuing 
with our simulation we get 



t 



j+3 



0'2, 0-3, f74, 1T5, (76, 0-7 = 21(/l + i) -5, 

as = Q{xi + 2), 

erg = 6(.T2 + 1), 

aw = 6{x3 + 1), 
(Til = 6, 



g21(/t+i)-5 



^ A, 



o'2, 0-3, (T4 = 21(/i + i) 
a5,ae,ar = 21{h + i) 
as = Q{xi + 3), 
erg = 6{X2 + 2), 
CTio = 6{X3 + 2), 
CTii = 6, 



7, 
7, 



^21(/.+.)-7/^2 
21(/i+i)-7 ^g21(ft+i-() + 10 



S, 



A, 



0-2,a'3,(T4 = 21(/i + i) 
(75, (76, 177 = 211 — 16, 

(78 =6(.Tl+4), 
(7g = 6(.T2 + 3), 

aio = 6{x3 + 3), 

(711 = 6, 



g21(/t+i)-8 



S, 



A. 



In configurations f^+s, and i^+s each of the simulated counters as, erg and crio are incremented. 
In configurations tj+e to tj+io the simulated counter as is decremented 4 times and the simulated 
counters ag and uio are each decremented twice. 



t 



j+6 



(72,(73 = 21(/l + i)-10, 

(74 = 21(/i + i) - 10, 

Cr5,0-6,cr7 = 21^ - 15, 

as = 6{xi + 4) + 3, 
0-9 = 6{x2 + 3) + 3, 
0-10 = 6(X3 + 3) + 3, 
(711 = 3, 



21(/i+j)-10 



^ S, 



21(/i+i)-10 yg21(/i+i-()+5 



■ s, 
A. 
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tj+7 : 

(72, (73 = 21(/l + Z) - 11, g21(h+i)-ll/g6 ^ 

(74, (75, (7(5, (77 = 21^ — 11, 

(78 = 6(a;i + 3) + 3, {s'^ys^s^ ^ s, 

ag = 6(X2 + 2) + 3, (s^)*.sV.s9 s, 

aio = 6(a;3 + 3) + 3, (sS)*.s^V*' ^ s, 

(711 =4, A, 



*j+8 : 



(72,(73 = 21(/i + i) - 13, 
(74, (75, (76, (77 = 21Z — 7, 

(78 = 6{xi + 2) + 2, 

(79 = 6{X2 + 1) + 2, 

o-io = 6(a;3 + 3) + 2, 
(711 = 3, 



21(/»+i)-13 1 g21{h+i-l)-G 



A, 



*j+9 • 

0'2,Cr3,C74,(75,(76,(77 = 21Z — 3, 

(78 = 6(a;i + 1) + 2, 

(79 = 6(X2 + 1) + 2, 

aio = 6(a;3 + 2) + 2, 

(7ii = 3, 



(s6)*s"/s2 ^ s, 
A, 



^■2, 0'3, (74, (75, (76, (77 = 21/, 

(78 = 6X1, 

(79 = 6(a;2 + 1), 

(710 = 6(0:3 + 1), 

(711 = 1, s ^ A, 

fTl2, CT13, = 21/1 - 29, {S^S'^/S^ S, 

CT15, cTie, fiy = 30. 

Note that at time tj+s that rule (s^)*s^'*/s^ — ?■ s will always be applicable as here X2 > (see the 
second line at the start of the proof). 
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The remainder of this siniulation is similar to the computation carried out at the end of the 
initiahsation process (sec the last paragraph of Section 4.0.6 and timesteps txi+x2+'i to ta;i+x2+i4/i.+3 
of the Section 4.0.5). Thus, after a further 14/i — 8 timesteps we get 

tj+14h+2 '■ 

(^2,o-3,a4,a5,a6,a7 = 21{h + + 

as = 6.T1, 

(Tg = 6(.T2 + 1), 

aio = 6{x3 + 1), 

<7l2, CIS, 0-14 = 21/l + 1, (s'^)*S^/s^ S. 

At timestep tj+i4h+2 the simulation of : DEC(l),qi,qk for Case 1 (a;i > 0) is complete. Note 
that a decrement on the value xi in counter ci is simulated by decreasing the value in cts from 
6{xi + 1) to 6x1. Note also that the encoding 21{h + I) + 1 of the next instruction qi has been 
established in neurons a2, as, 0-4, 0-5, ae and 0-7. Alternatively, if we have Case 2 (a;i = 0) then we 
get 

(72 =21(/i + i) -4, 

(73, 0-4, 0-5, 0-6, 0-7 = 21(/i + i) + 2, 



0-8 = 7, 

(79 = 6(a;2 + !) + !, 
(710 = 6(0:3 + !) + !, 

CTll = 1, 



(/)*sV/ 
(s'^)*sV/ 
s - 



A, 
' s, 

■ s, 
A, 



tj+2 : 



(72 =21(/i + i)-2, 

fs, 0-4, 0-5, (76, (77 = 21{h + i) +4, 

(79 = 62:2, 

o-io = 6X3, 

(711 = 1, 



21(h+i)-2 / 21(h+i-k)-l 



„21(ft,+5 



)+4 ^^21(ft,+i-fc)+5 



S, 



A, 



t 



■3+3 • 



0'2,C3,0'4,(75,(76,(77 = 21k, 
(78 = 6, 

(79 = 6(a;2 + 1), 
(710 = 6(a;3 + 1), 

(711 = 6, 

0'l2,(7i3,(7i4 = 21ft,- 8, 
O'lS, 0"16, CTl7 = 9- 



^ A, 



The remainder of this simulation is similar to the computation carried out at the end of the 
initialisation process (see the last paragraph of Section 4.0.6 and timesteps txi+a;2+4 to ixi+x2+i4/i+3 
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of the Section 4.0.5). Thus, after a further 14/i — 1 timesteps we get 

tj+14h+2 ■■ 

0-2, era, 0-4, 0-5, £76, (77 = 21(/l + fc) + 1, 

CTg = 6, 

erg = 6{X2 + 1), 
(710 = &{X3 + 1), 
a'12,cri3,0'l4 = 21/1 + 1. 

At time ^^+14^+2 the simulation of qi : DEC{l),qi,qk for Case 2 (xi — 0), is complete. Note that 
the encoding 21(/i + /c) + 1 of the next instruction q^ has been established in neurons a2, era, (74, 
(75, (76 and (77. 

4.0.9. ifaten^ 

If C3 enters the halt instruction qn at time tj then we get the following 

tr- 

0-2, era, C4, (^5 = 42/i + 1, 

(76,(77 = 42ft, + 1, 
(78 = 6(xi + 1), 

(79 = 6(a;2 + 1), 
(710 = 6(a;o + 1), 

tj+i : 

(72, (73, (74, (75 = 42/l + 1, 

(76, a7 = 42/1 + 2, 
£78 = 6(a;i + 1) + 4, 
£79 = 6(a;2 + 1) + 4, 
(710 = 6(a;o + 1) + 4, 
(711 = 4, 

ij+2 : 

£72, £73 =42/1 + 3, 
£74, £75 = 42/l + 3, 

C6 5 C7 = 42/i + 5, 
(78 = 6(xi + 2) + 2, 
(79 = 6(a;2 + 1) + 2, 
(710 = 6(a;o + 1) + 2, 
£711 = 5. 

Note that after time tj+2 we can ignore neurons (74, (75, (76 and (77 as there are no rules applicable in 
these neurons when the number of spikes is > 43/i + 3. The number of spikes in a2 and (73 does not 



„42/i+l 



/s S, 



g*^42/.+3/^ ^ 



(/)*s"/s2 ^ s, 
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decrease following timestep tj+2, and thus the rule s*s^'^^^^ /s — > s is applicable at each subsequent 
timestep regardless of the operation of neurons erg and cg. Thus, neurons cg and CTg may also be 
ignored as their operation has no effect on the remainder of the simulation. Note that in subsequent 
configurations we write (72, as > 42/i + 3 as there are more than 42/i + 3 spikes in each of these 
neurons. Thus we have 

(72, 0-3 ^ 42/i + 3, s*s'^'^^+^ls s, 

aio = 6xo + 2, {s'^ysys'^ ^ s, 

an = 8, 

tj+4 : 

(72, (73 ^ 42/l + 3, s*.s42h+3/_^ ^ 5^ 

aio = 6(xo-l) + 2, (s'^)*sV/ ^ s, 

CTii = 11, S^V*^ S> 

(72,(73 > 42/1 + 3, .s*.s''2''+3/s ^ s, 

(710 = 6(xo - 2) + 2, (s'^)*s7sS ^ s, 

(711 = 12. 

The rule (s^)*s*/s* — > s is applied in aio a further Xq — 2 times until we get 

tj+Xo+3 '■ 

(72, (73 > 42/i + 3, s*s^2fc+3^g _^ 

(710 =2, ^ A, 
an = S{xo - 2) + 12, 

a2, (73 ^ 42/1 + 3, s*s''2''+3/s _^ 

(710 =2, ^ A, 

(711 = 3(a;o - 2) + 14, {s^)*s^'^/s s. 

Recall from Section 2 that the output of an SN P system is the time interval between the first 
and second spikes that are sent out of the output neuron. Note from above that the output neuron 
(711 fires for the first time at timestep tj+4 and for the second time at timestep Thus, the 

output of is Xo the contents of the output counter C3 when C3 enters the halt instruction qh- 
If a^o = neuron an will fire only once. To see this, note that if = then s"^ X will be applied 
in neuron aio at time ij+3, and thus an will have 10 spikes (instead of 11) at time tj+4 and the 
rule s^'^ s will be; applied in an ending the computation. 

We have shown how to simulate arbitrary instructions of the form Qi : INC{l),qi and : 
DEC{l),qi,qk- Instructions that operate on counters C2 and C3 are simulated in a similar manner. 
Immediately following the simiilation of an instruction Ilj^^ is configured to begin simulation of the 
next instruction. Each instruction of C3 is simulated in 14/i + 2 timesteps. The pair of input values 
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origin ii(.'uiuiis 


t.ai 




iitnirolib 






0-2, 


0-3 


CT4,(T6,(76,(T7 


(T8, <T9,<Tl0, Cri2, <Tl3,<Tl4 


0-2, 


0-3, 


(T4 


as, ere, (77, (78 


(T9, <T10, cril 




"^2, 


0-8 


(T9, (Tio, cTll 






0-8, 


0-9 


CTIO, 0"11 






0"2, 


0-3 


<T4, iT5, fe, 




ClO 


0-2, 


0-3 


CT4,(T6,(76,(T7 


CTll 


0'12,0'13, fl4 


CIS 


, CIS, fir 




flSi^ie; cri7 


0-2, 


0-3 


CT4,(T6,Cr6,C7 


0'12,<Tl3,Cri4 



Table 3 

This tabic gives the set of synapses of the SN P system n|-,^ . Each origin neuron (Ji and target neuron aj that appear 
on the same row have a synapse going from ffj to aj. 




—S 



Fig. 4. Finite state ma<;hine G decides if there is any rule applicable in a neuron given the number of spikes in the 
neuron at a given time in the computation. Each s represents a spike in the neuron. Machine G' keeps track of the 
movement of spikes into and out of the neuron and decides whither or not a particular rule is applicable at each 
timestep in the computation. +s represents a single spike entering the neuron and — s represents a single spike exiting 
the neuron. 

{xi,X2) is read into the system in xi + X2 + 14/i + 3 timesteps and sending the output value Xo 
out of the system takes + 4 timesteps. Thus, if C3 completes it computation in time t then 11^^ 
simulates the computation of C3 in linear time 0{ht + xi + X2 + Xo)- □ 

5. Lower bounds for small universal SN P systems 

In this section we show that there exists no universal SN P system with only 3 neurons even 
when we allow the input technique to be generalised. This is achieved in Theorem 4 by showing 
that these systems are simulated by log-space bounded Turing machines. Following this, we show 
that if we generalise the output technique we can give a universal SN P system with extended rules 
that has only 3 neurons. As a corollary of our proof of Theorem 4, we find that a universal SN P 
system with extended rules and generalised input and output is not possible with 2 neurons. 
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In this and other work [15,18] on small SN P systems the input neuron only receives a constant 
number of spikes from the environment and the output neuron fires no more than a constant number 
of times. Hence, we call input standard if the input neuron receives no more than y spikes from 
the environment, where y is a constant independent of the input (i.e. the number of Is in its input 
sequence is < y). Similarly, we call the output standard if the output neuron fires no more than x 
times, where a; is a constant independent of the input. Here we say an SN P system has generalised 
input if the input neuron is permitted to receive ^ n spikes from the environment where n G N is 
the length of its input sequence. 

Theorem 4 Let H he any extended SN P system, with only 3 neurons, generalised input and stan- 
dard output. Then there is a non- deterministic Turing machine Tn that simulates the computation 
o/n in space O(logn) where n is the length of the input to 11. 

PROOF. Let n be any extended SN P system with generalised input, standard output, and 
neurons (Ti, (T2 and 173. Also, let x be the maximum number of times the output neuron 173 is 
permitted to fire and let q and r be the maximum value for b and p respectively, for all E/s^ — > 
sP; d in H. 

We begin by explaining how the activity of (T3 may be simulated using only the states of Tn (i-e. 
no workspace is required to simulate 0-3). Recall that the applicability of each rule is determined by 
a regular expression over a unary alphabet. We can give a single regular expression R that is the 
union of all the regular expressions for the firing rules of (T3. This regular expression R determines 
whither or not there is any applicable rule in as at each timestep. Figure 4 gives the deterministic 
finite automata G that accepts L{R) the language generated by R. During a computation we may 
use G to decide which rules are applicable in Ui by passing an s to G each time a spike enters 0-3. 
However, G may not give the correct result if spikes leave the neuron as it does not record spikes 
leaving o-j. Thus, using G we may construct a second machine G' such that G' records the movement 
of spikes going into and out of the neuron. G' is construct as follows: G' has all the same states 
(including accept states) and transitions as G along with an extra set of transitions that record 
spikes leaving the neuron. This extra set of transitions are given as follows: for each transition on 
s from a state gi to a state gj in G there is a new transition on — s going from state gj to gi in G' 
that records the removal of a spike from 0-3. By recording the dynamic movement of spikes, G" is 
able to decide which rules are applicable in as at each timestep during the computation. G' is also 
given in Figure 4. To simulate the operation of as we emulate the operation of G' in the states of 
Tn- Note that there is a single non-deterministic choice to be made in G' . This choice is at state 
gu if a spike is being removed (— s). It would seem that in order to make the correct choice in this 
situation we need to know the exact number of spikes in as- However, we need only store at most 
u + yq spikes. The reason for this is that if there are ^ u + yq spikes in as, then G' will not enter 
state gu-i again. To see this, note that as spikes a maximum of y times using at most q spikes 
each time, and so once there are > u + yq spikes the number of spikes in as will be > u — 1 for the 
remainder of the computation. Thus, Tn simulates the activity of as by simulating the operation 
of G' and encoding at most u + yq spikes in its states. 

In this paragraph we explain the operation of Tn. Following this, we give an analysis of the space 
complexity of Tq. Tq has 4 tapes including an output tape, which is initially blank, and a read 
only input tape. The tape head on both the input and output tapes is permitted to only move 
right. Each of the remaining tapes, tapes 1 and 2 simulate the activity of the neurons <Ji and 1T2, 
respectively. These tapes record the number of spikes in cti and a2- A timestep of H is simulated 
as follows: Tn scans tapes 1 and 2 to determine if there are any applicable rules in ai and a2 at 



28 



that timestep. The applicability of each neural rule in 11 is determined by a regular expression 
and so a decider for each rule is easily implemented in the states of Tu- Recall from the previous 
paragraph that the applicability of the rules in CT3 is already recorded in the states of Tu- Also, Tn 
is non-deterministic and so if more than one rule is applicable in a neuron Tn simply chooses the 
rule to simulate in the same manner as 11. Once Tn has determined which rules are applicable in 
each of the three neurons at that timestep it changes the encodings on tapes 1 and 2 to simulate 
the change in the number of spikes in neurons cri and (T2 during that timestep. As mentioned in the 
previous paragraph any change in the number of spikes in (T3 is recorded in the states of Tn. The 
input sequence of 11 may be given as binary input to Tn by placing it on its input tape. Also, if 
at a given timestep a 1 is read on the input tape then Tn simulates a spike entering the simulated 
input nciiron. At each simulated timestep, if the oiitpiit neuron (J3 spikes then a 1 is place on the 
output tape, and if CT3 does not spike a is placed on the output tape. Thus the output of If is 
encoded on the output tape when the simulation ends. 

In a two neuron system each neuron has at most one OTit-going synapse and so the number of 
spikes in the system does not increase over time. Thus, the total number of spikes in neurons ai 
and a2 can only increase when (J3 fires or a spike is sent into the system from the environment. 
The input is of length n, and so cti and (T2 receive a maximum of n spikes from the environment. 
Neuron (T3 fires a total of y times sending at most r spikes each time and so the maximum number 
of spikes in ai and CT2 during the computation is n + 2ry. Using a binary encoding tapes 1 and 2 
of Tn encode the number of spikes in ci and <T2 using space of log2(n + 2ry). As mentioned earlier 
no space is used to simulate 0-3, and thus Tn simulates 11 using space of O(logn). □ 

It is interesting to note that with a slight generalisation on the system in Theorem 4 we obtain 
universality. If we remove the restriction that allows the output neuron to fire only a constant 
number of times then we may construct a universal SN P system with extended rules and only 
three neurons. Here we define the output of an extended SN P system with generalised output to 
the time interval between the first and second timesteps where exactly x spikes are sent out of the 
output neuron. 

Theorem 5 Let C2 be a universal counter machine with 2 counters that completes it computation 

in time t to give the output value Xq when given the input value xi . Then there is a universal extended 
SN P system 11^^ with standard input and generalised output that simulates the computation of C2 
in time 0{t + xi + Xo) and has only 3 neurons. 

PROOF. A graph of 11^^ is constructed by removing the output neuron (75 from the system Hc^ 
given in the proof of Theorem 2 and making 0-3 the new output neuron of 11^^ . The rules for 11^^ 
are given by the first 3 rows of Table 5 and a diagram of the system is obtained by removing neurons 
(74 and (75 from Figure 1 and adding a synapse to the environment from the new output neuron 
0-3 . The operation of 11^^ is identical to the operation of Tic^ with the exception of the new output 
technique. The output of 11^^ is the time interval between the first and second timesteps where 
exactly 2 spikes are sent out of the output neuron 0-3. □ 

Prom the third paragraph of the proof of Theorem 4 we get the following immediate corollary. 

Corollary 1 Let 11 he any extended SN P system with only 2 neurons and generalised input and 
output. Then there is a non- deterministic Turing machine Tn that simulates the computation of 11 
in space O(logn) where n is the length of the input to 11. 
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6. Conclusion 



The dramatic improvement on the size of earher smaU universal SN P system given by Theorems 1 
and 3 is in part due to the method we use to encode the instructions of the counter machines our 
systems simulate. In the systems of Paun and Paun [15] each counter machine instruction was 
encoded by a unique set of neurons. Thus the size of the system is dependant on the number of 
instructions in the counter machine being simulated. Some improvement was made by Zhang et 
al. [18] by showing that certain types of instructions may be grouped together. However, the number 
of neurons used by the system remained dependant on the number of instructions in the counter 
machine being simulated. In our systems each unique counter machine instruction is encoded as 
a unique number of spikes and thus the size of our SN P systems arc independent of the number 
of instruction used by the counter machine they are simulating. The technique of encoding the 
instructions as spikes was first used to construct small universal SN P systems in [14]. 

The results from Theorems 2 and 4 give tight upper and lower bounds on the size of the smallest 
universal SN P system with extended rules. Thus in Theorem 2 we have given the smallest possible 
universal SN P system with extended rules. The results from Theorem 5 and Corollary 1 give tight 
upper and lower bounds on the size of the smallest universal SN P systems with extended rules 
and generalised output. Thus, Theorem 5 gives the smallest possible universal SN P system with 
extended rules and generalised output. 

The lower bounds given in Theorem 4 are also applicable to standard SN P systems and thus give 
a lower bound of 4 neurons for the smallest possible standard system that is universal. However, 
when compared with extended systems the rules used in standard SN P systems are quite limited, 
and so it seems likely that this lower bound of 4 neurons can be increased. Note that here and 
in [15,18] the size of a universal SN P system is measured by the number of neurons in the system. 
However, the size of an SN P system could also be measured by the number of neural rules in the 
system. 
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neuron 


rules 


cri 


gSh+l IgSh ^ gSh^ g8h+2 l^Sh+l _^ g8h+l _ g6h+l _^ g4h+4^ -j. A, S A 




gl6h+4i _^ gl2h+4l _ glOh+4i ^ g4(h+l) ^ if Kh 




gl6h+4i ^ ^12^+5^ glOh+4i _^ g4h+5^ ii I = h 




gSh+4i ^ g4(h+k) ^ iffcT^/l 




gSh+4i _^ g4h+5^ iik = h 




^gShygSh+l/gSh _^ (^g8hygSh+2/g8h+2 _^ g2h 




^g8h-^*g4(h+i)/g4{h+i) ^ g4(h+i) if g. . /jvc(i) e {Q} 




(^g8h-j*g4{h+i)/g8h+4(h+i) _^ g6h if g. . /ivC(a;) € {Q}, x ^ 1 




(^8h)*^16h+4(h+i) /^12h+4i ^ ^6h+4i jf ^_ . DEC(l) G {Q} 




g8h+4(h+i) /g4(h+i) _^ s4,(h+i) if . /)£'C'(1) g {Q} 




(^gShy g4(h+i) /g4ih+i) _^ g2h if g. . DEC{x) e {Q}, X 7^ 1 


0-3 


s'^/s^S, sS'^+l/s^'' (s^hyg8h+3/g8h+3 ^ g2h^ (^8h)* ^20h+5/gm _^ ^2 




(58h)*5l6h+5/58/» _^ 5^ 58^+5 _^ ^2^ ^12h+5 _^ ^2 




(S^hyg4(h+i)/g4(h+i) _^ g4(h+i) if g. . /JVC(2) e {Q} 




(^gShy g4(h+i) ig8h+4(h+i) _^ g&h if . J^C{x) e {Q}, X + 2 




(^8h).^16h+4(h+i) jg\2h^4i _^ g(,h+4i if . DEC{2) e {Q} 




^8?.+4(h+i) /34(/i+i) _^ 54(/i+i) if . DEC(2) g {Q} 




(^8h).^4(h+i)/^4(h+i) _^ ^2^ if g. . DECix) G {Q}, x^2 


0-4 


g8h+l/g8h _^ g8h-l^ s^^+^/s^'' ^ S^h 




(g8h)*g4(h+i)/g4(fe+i) _^ 54(/i+i) if ^. . /JVC(3) 6 {Q} 




(^8h).^4(h+i) ^^8h+4(h+i) _^ g6h if . JJ^C{x) e {Q}, X + Z 




(^8h).^16h+4(h+i) /^12h+4i _^ ^6h+4i jf ^. . f)EC{3) € {Q} 




^8J.+4(h+«)/54(h+i) ^ 34(h+i) if . DECiS) e {Q} 




(^8h).^4(h+i)/54(?,+i) _^ ^2ft if ^. . DECix) e {Q}, a: ^ 3 




s -> A, s^^ -> A, A, s^C'+i) ^ A, s6/>+4i _^ -> s 



Table 4 

This table gives the rules for each of the neurons of flc^ . 
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neuron 


rules 


cri 


g8h+l^g8h ^ gSh^ g8h+2^g8h-l _^ s4/i+3^ g8h+3 ^ 


A, s — > A 




gl6/i+4i _^ gl2h+4l ^lQh+4i _^ ^4(/i+i) ii I <C h 






lOd^c; inii^/i^' .n, , 
gl6fi+4t _j. ^12h + 5^ glOft+4i _^ ^4h+5^ if i = ft, 






58/1+4* _j. ^4(h+fc)^ iffe^/l 






38/»+4j _j. 34/1+5^ ififc = ft 






(^g8hyg4(h+i)/g4(h+i) _^ 54(/,+i) if q. . INC(1) e {Q} 






(g8h)*54(h+i)/58h+4(/i+i) _^ ^12/1 if . 7JVC(2) e {Q} 






(58h)*5l6h+4(h+i) /5l2h+4i _^ g6/i+4j if . DEC{1) 6 {Q} 






^8/1+4(^+1) /54(h+») _^ g4(h+i) if . £)EC(1) e {Q} 






(s8'')*s4(h+i) /s4{h+i) ^ ^4(1 if ^. . DEC {2) e {Q} 




0-3 








(s8h).5l6/»+6/58/» _^ 5^ 58^+5 _^ ^m+B _^ ^2 






(^S^hy g4(h+i) 1 g4(h+i) ^ g4{h+i) if . 7JVC(2) 6 {Q} 






(58h).54{h+i) Ig8h+4{h+i) _^ gl2h if ^. . 7JVC(1) e {Q} 






(58h).5l6h+4(h+») |gl2h+4^ _^ ^bh+4i if . UEC(2) S {Q} 






58h+4(h+»)/54(h+») ^ 54(h+«) if g. . DEC {2) G {Q} 






(58h)*54{h+i) /34{/i+i) ^ ^4/1 if ^, . DEC{1) G {Q} 




<T5 


s*'' A, si2h. _^ s A, s*'' A, s6h+4j _^ 





Table 5 

This table gives the rules for each of the neurons of IIc^ • 
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neuron 


rules 








s41/s-).s, S^V^^^-^S, s22/s3^S, 






321(h+i 


-i/s^^S, s^lC'+O-S/gS g21(h+i)-7 /g2 _ 


> 5, g2X(h^i)-%lgZ 5 




g2\(h + i 


~W/g5 ^ 3_ g42h+l Jg ^ g^ g,gi2h+3lg _^ g 






g1\{h + i 


+1/5" s if 9i : INC G {Q} 






g21(h+i 


-2lg2l(h+i-l)+l ^g if . INC{1) e {Q} 






g21ih+i 


-2/^2i(h.+i-o+2 _^ 5 if . /jvC(x) e {Q}, x ^ 1 






g2l(h+i 


-2/^21(a+i-fc)-l _^ g if ^. . g 






g21(h+i 


+ 1/5-°^ if (7i : D£;c(i) e {Q} 






g21(h+i 


+Vs'5 ^ s if ; DEC{x) e {Q}, X 7^ 1 






g21(h+i 


-11/36^^ a q^ : DEC(l) € {Q} 






g21{h+i 


-"/■s^-^s if gi : DEC(a;) G {Q}, x 7^ 1 






g21{h+i 


-i3/52i(a+i-0-6 _^ 5 if g. . DEC(1) e {Q} 






g21ih+i 


-13/^21{h.+<-0-7 ^g if . g {Q}, X 7^ 1 






s-ll/s - 




sSlC'+O-l/sB ^ g 




g21(h+i 


s21('^+')-5/s^ ^'^> ^21(h+z)-7/^2 


s, g2Hh+i)-8/g3 ^ g 




g21(h+i 


-11/36 5 j,21{/i+i)-13/^21(h+i-0-6 _^ g21{h+i)+4/^21(h+i-fc)+5 _^ ^ 




g42h+l. 








g21(h+i 


+i^g2i(h+i-i)+6 _^g if . /ivc(i) e {Q} 






g21(h + i 


+ i/s^ ^ A if qi : INC{x) 6 {Q}, x ^ 1 






g21{h + i 


-2/^,21(fe + i-0+2 ^ ^ if g. . J]^C{x) G {Q} 






g21(h + i 


-2/^,2i(fe+i-fc)-l ^ ^ ii ; DEC & {Q} 






g21{h+i 


+ ^- s if 9i : DEC{x) G {Q}, a; 7^ 1 






g21(h+i 


-10/55 _^g if g. . DEC(l) e {Q} 






g21(h+i 


-10/g21(h+i-l)+5 ^g if g. . G {Q}, X + 1 





Table 6 

This table gives the rules for neurons ui to 1T3 of 



34 



neuron 



rules 



S^^/s^S, s^^/s^^S, s44/^31_^^_ 516/^3^^^ ^21(h+i)-l/gS 

s2^h + ^) + 5/^ll _^ ^21(fe + »)-5/^3 ^ ^2^h + ^) -8 / ^3 ^ ^ 

^21(h+i)-10 ^^21(h+i-0+5 _^ ^21(h+j)+4 ^^21()i+i-fc)+5 _^ ^42h+l _^ ^ 

g21{h+i)+l^g21(h+i-l)+6 if q. . I]^fC{x) £ {Q}, X^3 

g2i{h+i)+i/g4 if ^. . iNC{3) e {Q} 

g2iih+i)-2/g2i(h+i-i)+2 if . ir^c{x) e {Q} 

^21(h + i)-2 /^21(h + ^-k)^l if . g ^gy 

g21(h+^) + l/g6 if g, . DEC{3) e {Q} 

g2Hh+,)-7/g2i(h+i-l)+io ^ g if q. . DEC{3) e {Q} 
g2i(h+i)-7/g2 if . oECCa;) G {Q}, x 7^ 3 

^21(ft+i)+6/gll _^ ^21(a+i)-5/53 _^ 5^ g21(h+i)-7/g21(h+i-l)+10 _^ ^ 

g21(h+i)+4 yg21(;i+i-fe)+5 _^ ^42h+l ^ s 

^,2i(fe+i)+i/^2i(?i+«-o+6 if g_ j^c* e {Q} 

^21(h+i)+5/5ll g21(h+i)-5/g3 ^ g^ g21(h+i)-7/g21(h+i-l)+W ^ g 

g21(h+i)+4 1 g21(h+i-k)+5 _^ g 

g21(h+i)+l ig2l(h+i-l)+6 _^g if . g {Q} 



0-8 


(s6)*s"/s« ^ S, {s6)*sl3/s ^ S^^A, 
(s'O'sVs*^ ^ s. (s6)*sVs-' ^ s, 


0-9 


(s6)*sio/se s, ^- s, (s6)*si4/52 _^ 5, 

58 ^- A, (s6)*s9/s9 s, 


flO 


(s6)*si0/s6 s, (s6)*s"/s® -> s, (s«)*s''/s'' -> S, 




s^'^-A, s^^X, s-i-A, siV«->s. (s3)*si4/s s, 

s4 ^ A, s2 ^ A, s3 ^ ^10 _^ 5 


0-12, cri5 


(s'^)*s*/s'^ s, s s, 


"■13, 014, (716, cri7 


(s3)*s4/53 _^ s ^ A, 



Table 7 

This table gives the rules for neurons 1T4 to 1T17 of nj- 
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